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An explicit solution of time-homogeneous pure birth branching processes is 
described. It gives alternative extensions for the negative binomial distribu- 
tion (branching processes with immigration) and for the Furry- Yule distribution 
(branching processes without immigration). 



1 Introduction 

The Furry- Yule (FY) and negative binomial (NB) distributions are widely used 
in phenomenological studies of multiplicity distributions at high energies. Re- 
cently, several generalizations of the negative binomial distribution based on 
the perturbaliye quantum chromodynamics have appeared in the literature (s££ 
for example Era). The aim of this note is to present another generalization uu 
of the Nliand-FY distributions based on simple time- homogeneous branching 
processestfQ'Q'cl more simple than the branching processes used in QCD. The 
NB and FY distributions occur in time-homogeneous branching processes with 
allowed transition 1 — > 2 and one can extend them allowing higher order tran- 
sitions 1 — > n with n > 2. It may happen that these processes can be useful 
in the description of the hadronization stages still not fully understood in the 
framework of the QCD. 

In section 2, we describe the distribution for the pure birth branching 
process with allowed higher-order transitions that can be used as an extension 
of the FY distribution. In section 3, the distribution for the branching process 
with immigration □ is presented. A discussion is given in the last section. 



2 Distributions for general pure birth branching processes 

A branching process with continuous evolution parameter t is determined by 
the rates a n for the transition ( "splitting" ) of one particle into n particles with 
all particles subsequently evolving independently. For a pure birth branching 
process ao = 0. The FY distribution occurs when only the oci is not equal to 
zero. The probability generating function m(x, t) = Y^=oPri{t)x n satisfies the 
forward Kolmogorov equation&0 ^ = f(x)^ with f(x) — X)^2 a nX n — ax 
and a = a n- The Taylor expansion of this equation leads to the following 
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system of equations for the probabilities p n 

dpi m 

"dt" = ' (1) 

and for n > 1 

, n— 1 

dp™ • ^ 

= 2^^n-j+iPj - nap n . (2) 

For the case of N initial particles with Pn (0) — <5jv„, the solution of this 
system has the following formQ 

n 

PI N) = E , (3) 

with the coefficients 7rj„ obeying the recursion: 

n—j 

(-iH£° = E(«- Ofc-iS-o ■ (4) 

Here, pi = exp(— at) and 6/ = ai+i/a. This recursion starts from 7rjy^ = 1 
and the coefficient 7Tnn can be found from the relation tt^P = — X^J^jv 71 jn ■ 
One can calculate the coefficients 7rj„ using the concept of the Koenigs 
function (seei'0 and the references therein). For the branching process starting 
from N particles, the Koenigs function is defined as the limit: 

= = n fim = E = E ^ ■ (5) 

The recursion (|]) leads to the following recurrence for the coefficients Kj , N = 
1 J A • 1..V • 2.... : 

j-AT 

(i-iv)^ = E(i-0^^) ■ (6) 

It is convenient to denote «4+ n =_f„(a;), then to(x) = 1 and the polynomials 
t n (x) are defined by the recursionQ 

n 

nt n (x) = )Xx + n- l)h t n -i{x) . (7) 
i=i 
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The /tj. in terms of t n (x) is equal to tj-N(N). 

The remarkable property of the Koenigs function is that it satisfies the 
functional Schroder equation K{m) — p\K(x). This can be easily obtained 
from the forward Kolmogorov equation dm /At = f(m). Its integration gives 
also the expression 

K(x) =exp(-a f 4r^\ • (8) 



/(«) 

It is convenient to introduce the function Q{x), the inverse of the Koenigs 
function. Then, the inversion of the Schroder equation gives the functional re- 
lations m(x,t) — Q(p\K(x)) and m (x,t) = Q N {p\K{x)). The coefficients of 
the Taylor expansion for Q N (x) = Qj can be found using the following 
relation: = ^K ~Jf ■ One can show that this relation is equivalent to the 

Biirmann-Lagrange series for the inverse function (see Appendix B in0 and 
the references therein). In terms of t n (x) it gives: Qj — yij-jv(— j). Finally, 
comparison of equation ([?]) with the Taylor expansion in x of the Q N (p\{K (x)) 
leads to the following expression! for the coefficients in terms of t n (x): 

t$° = Qf = K t H)t {j) . (9 ) 

In some cases, the functions K(x) and Q(x) can be found explicitly, for 
example: 
for bi = 1 



K{x) = -^- , Q(*) = -J- ; (10) 
1 — x 1 + x 



for 6iv = 1 



K{ - x) = (i - x N y/ N ' = (i + x N y/ N ; 

for bi + b 2 = 1 

K{x) = — , Q(x)=? . (12) 

(1 - x) 1+b 2 (1 + b 2 x) r+^i 

In the general case, the expressions for Q(:c) and K(x) are quite complicated. 

The solution for non-critical branching processes with non-zero ao is an 
infinite series in exp(— a't), see □. 
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3 Solution for branching processes with immigration 

For the branching processes with immigration, there is an additional external 
source of particles appearing generally in clusters of j particles with the dif- 
ferential rates Pj ( Pj = b ). The probability generatingJjunction for the 
process starting with zero particles at t — can be written EttI3 as 

M(x,t) = exp^ g(m(x,T))dT\ (13) 

with 

oo 

g{x)=Y J ^ i -b , (14) 
where m(x,r) is the solution for the underlying branching process without 



immigration. For the underlying pure birth branching process. Equ. (13) 
leads to the following expressionQ 

M{x,t) = exp (-fa) exp C n (t)xA (15) 



with 



Let us denote 



a=l j=i 



exp (J2 C n x n ) =l + Y d V n x n . (17) 

n— 1 n— 1 

Then, the final probability P n (t) = exp (— bt)V n (t). The coefficients V n (t) can 
be calculated using the recursive relation: 

n 

nV n = Y,jC j V n - j . (18) 

7=1 

iaJmnffipJn combinatorics and is used, for example, in the study 
iBHota The coefficients C n (t) are known as the combinants of 
the multiplicity distribution P n (t) . It is of interest to note the formal analogy 
between Equ. ( |l7j ) and the exponential representation for the K(x) (JsJ) . 

The branching process with immigration leads to the NB distribution when 
only Pxi cx-2 and czq are not equal to zero. The resulting distribution is the 
Poissonian when only /3% is non-zero. 



This relation 
of combinants 
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4 Discussion 



In this note we have described explicit expressions for the probability distribu- 
tions in the branching processes with the parameters equal to the differential 
rates oti+\ and /3j. These distributions can serve as an extension for the NB 
and/or FY distributions when at least one of the coefficients with i > 1 is 
non-zero. It seems natural to assume that these parameters have no energy 
dependence and therefore the energy dependence is concentrated in one param- 
eter pi(t) = cxp(— at). It is of interest to note that the exponential dependence 
on t corresponds to the power-law dependence on the mean multiplicity (n), 
since (n) has linear dependence on exp(/'(l)i). 

We have applied the distribution with non-zero «2 and to the e + e _ 
multiplicity data in the paperd It has been assumed that a fixed number N 
of sources of particle production is produced at some initial stage of the in- 
teraction. These sources develop independently of each other according to the 
time-homogeneous branching process, producing intermediate neutral clusters. 
Finally, these clusters decay into a pair of charged or neutral hadrons with 
probabilities e or (1 — e), respectively, as in the Goulianos modelO The data 
at cm. energies yfs above 20 GeV are well described with TV = 7, with energy 
dependent parameter p\ and with the parameter e fluctuating near 0.677. It 
has been shown that the parameter p\ has power-law dependence on s/s, as 
noted earlierta It has been also shown that 03 is consistent with zero, but val- 
ues of the ratio a^/a of the order ~ 0.05 are not excluded at the present level 
of statistics. It is of interest to note that the values of e near 2/3 are consistent 
with zero isospin of the intermediate-neutral clusters and are different from the 
value e = 0.5 used by K. Goulianos til 

It is a challenging task to derive explicit expressions for the polynomials 
t n (x) from the recurrence (j?]). This can be done easily when only one of the 
hi is not equal to zero, for example, for b\ = 1 the t n (x) = x ( x + 1 )---( x + n - 1 ) _ 

It is of interest to look for the singularities in the complex plane of the 
probability generating functions for the time- homogeneous branching processes 
(as pursued by I. Dremint3 and E. Dc WolfE3). These singularities are con- 
nected with the singularities of the K(x) and Q{x) and, therefore, with the 
roots of the equation f(x) — 0. For example, for the process with = 1, the 
poles are: 

- - eXP( ^» (19) 



(1-0"" 



(n = 1, 2, N). They form a regular polygon and move with t to the circum- 
ference with unit radius. 



5 



Acknowledgments 



I am grateful to S. Hegyi for the opportunity to give this talk and I would like to 
thank the organizers of the Symposium for providing a beautiful surrounding. 

References 

1. S. Hegyi, Phys.Lett. B 387, 642 (1996); ibid. B 388, 837 (1996). 

2. R. Botet and M. Ploszajczak in Proc. 28th Int. Symp. on Multi-particle 
Dynamics, Delphi, eds. Antoniou et al. (World Scientific, Singapore, 
2000) p.217. 

3. O.G. Tchikilev, Phys.Lett. B 471, 400 (2000); erratum, Phys.Lett. B 
478, 459 (2000). 

4. O.G. Tchikilev, Phys.Lett. B 491, 36 (2000). 

5. For a recent review see: R.C. Hwa in Hadronic Multiparticle Production, 
ed. by P. Carruthers (World Scientific, Singapore, 1988) p. 556. 

6. N.A. Dmitriev and A.N. Kolmogorov, Dokl. Acad. Nauk SSSR 56, 7 
(1947). 

7. B.A. Sevastyanov, Uspekhi Math. Nauk 6, 47 (1951). 

8. M.S. Bartlctt in An introduction to stochastic processes with special ref- 
erence to methods and applications, Cambridge Univ. Press, Cambridge, 
1955. 

9. M. Kuczma in Functional equations in a single variable, PWN - Polish 
Scientific Publishers, Warszawa, 1968. 

10. M. Kuczma, B. Choczewski and R. Ger in Iterative functional equations, 
Cambridge Univ. Press, 1989. 

11. J.F. Traub in Iterative methods for the solution of equations, Chelsea 
Publishing Company, New York, 1982. 

12. P.V. Chliapnikov and O.G. Tchikilev, Phys. Lett. B 235, 347 (1990). 

13. M. Gyulassy and S.K. Kauffmann, Phys. Rev. Lett. 40, 298 (1978). 

14. S.K. Kauffmann and M. Gyulassy, J. Phys. A 11, 1715 (1978). 

15. A.B. Balantckin and J.E. Soger, Phys. Lett. B 266, 231 (1991). 

16. S. Hegyi, Phys. Lett. B 309, 443 (1993); ibid. B 318, 642 (1993); ibid. 
B 463, 126 (1999). 

17. K. Goulianos, Phys. Lett. B 193, 151 (1987). 

18. M. Biyajima, T. Osada and K. Takei, Phys. Rev D 58, 037305 (1998). 

19. I.M. Dremin, JETP Letters 60, 764 (1994), Pis'ma Zh. Eks. Teor. Fiz. 
60, 757 (1994). 

20. E.A. De Wolf in Proc. of the XXIV Int. Symp. on Multiparticle Dy- 
namics, Vietri sul Mare, Italy, 1994, eds A. Giovannini, S. Lupia and 
R. Ugoccioni (World Scientific, Singapore, 1995) p. 15. 



6 



